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Abstract 

Fermion production by the electric field of a charge on de Sitter expanding universe is analyzed. 
The amplitude and probability of pair production are computed. We obtain from our calculations 
that the modulus of the momentum is no longer conserved and that there are probabilities for 
production processes where the helicity is no longer conserved. The rate of pair production in an 
electric field is found to be important in the early universe when the expansion factor was large 
comparatively with the particle mass. 

PACS numbers: 04.62,+v 



1 



I. INTRODUCTION 



Pair production by electric field is known as the Schwinger effect 25(, and even if this 
phenomenon was not yet observed, it might be possible to be experimentally proven in the 
next years. The method employed in obtaining this result is based on the non-perturbative 
structure of the QED vacuum physics, when strong fields are considered. However, pair 
production by a weak electric field is in principle possible if this field is coupled with strong 
gravity. This kind of conditions were achieved only in the early universe. In the present 
paper we want to prove that pair production by the electric field of a charge could have 
arisen in the early universe due to the large values of the expansion factor. 

The problem of fermion production by a constant electric field on an expanding universe 
was the subject of a few investigations 19, [2^, [^j . The main results of the mentioned 



papers show that there are indeed nonvanishing rates for pair production in an electric field. 
However, pair production in the early universe was not addressed in 19, Q, l^] due to some 



constrains which are imposed to the particle mass, which must be larger than the expansion 
factor of the space. 

It is important to specify that in our calculations we will use the electric field produced 
by a point charge (or a distribution of charges) in the de Sitter universe, which is a time 
dependent quantity. To the best of our knowledge, this kind of electric field was not consid- 
ered in the investigations related to the pair production in an expanding de Sitter universe. 
The work of L. Parker js lo| shows that we must expect a considerable amount of pair 
production from vacuum in the early universe. Then we must expect the phenomenon of 
pair production to be amplified in the presence of an external electric field coupled with 
gravity. 

In the present paper, we want to prove that the electric field produced by a charge, 
or a distribution of charges can produce large amounts of fermions in the early expansion 
conditions. Our method is based on perturbations and was used in a number of investigations 



for calculating QED processes in de Sitter universe [16 



181 ] . To the best of our knowledge, the 



method based on perturbations was not used to calculate the probability of pair production 
in electric fields on de Sitter space. We will compute the amplitude corresponding to the 
process vacuum — > e~ + e + in the presence of an electric field. The time reversed process can 
be also calculated, and represents the process in which the pair is absorbed by the electric 
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field. As we know, a point charge can produce an electric field with intensity decreasing 
with the square of distance to the source (charge). In flat space, this kind of electric field 
doesn't produce pairs of fermions. However, in the early universe, this phenomenon can not 
be excluded. 

The second section begins with a short review of the free fields theory. In the third 
section we compute the amplitude and probability of pair production in an electric field and 
we analyze the physical consequences of our calculations. The fourth section is dedicated to 
the calculation of the total number of produced fermions and in section five we present our 
conclusions. 



II. PRELIMINARIES 

The theory of the free Dirac field on de Sitter space is well known and in the Appendix 
is given only the form of the fundamental solutions, which have a definite momentum and 
helicity. These solutions will help us calculate the form of our amplitude. 

We will focus in this section in establishing the form of the electric field in the de Sitter 
metric jjj : 

ds 2 = dt 2 - e 2ult dx 2 = —^—;{dt 2 c -dx 2 ), (1) 

[L)t c ) 

where u is the expansion factor (u > 0) and the conformal time is defined as t c = — • The 
form of the line element allows us to choose the simple Cartesian gauge with the nonvanishing 
tetrad components 

4 = e"*"; el = 8Ke- ut ; e° = ; e) = 5}e ut (2) 

so that = e v ~e v and have the orthonormalization properties e^ep = 77^7, 
^"ep = Sp, with respect to the Minkowski metric r\ = diag(l, —1, —1, —1). 

The theory of the free electromagnetic field on de Sitter space is simple, if we recall that 
the de Sitter metric is conformal with the Minkowski metric. At a conformal transformation 
the metric tensor g a p can always be expressed in the form g a p = Q r) Q /3 , g al3 = Q~ l r] al3 , where 
r] a /3 is the Minkowski metric and Q = ,^ ^ . Using this transformation, it was shown that 
the four vector potential in de Sitter space can be expressed in terms of the corresponding 



four vector potential on Minkowski space as follows 
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A,s = A, M = (3) 
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The electric field of a charge will be under the expansion conditions a time-dependent quan- 
tity. This can be seen from the expression of the electromagnetic field tensor in the de Sitter 
metric. The electromagnetic field tensor in de Sitter space can be expressed in terms of 
the electromagnetic field tensor from Minkowski space, in order to preserve the conformal 
invar iance: 

F^uS = FfivM ,Fg U = Vl 2 F£f, (4) 

Then in a local Minkowski frame, equation fl4]) becomes, Fg V = e^e^F^ = e^e^f2~ 2 F%f . 
Now, if we consider the electric field produced by a charge Q, or a distribution of charges 
Qii we w iH find in de Sitter space, that at distance \x\ from the charge, the intensity of 
the electric field is: 

where n is the unit vector which gives the direction and orientation of the intensity of electric 
field. It is now simple to establish the form of the vector potential that will be associated to 
this electric field: A = 2 J^ 2 e~ 2wt n. Then, in a local frame the potential vector, = e^A v 
will be: 

1 = -JL-e-^fi. (6) 



2oj\x 



This expression will be further used to calculate the amplitude of pair production by an 
electric field in de Sitter universe. 

III. TRANSITION AMPLITUDE AND PROBABILITY 

In this section we will calculate the transition amplitude and probability, for pair pro- 
duction in an electric field and we will explore the physical consequences of our result. The 
form of the transition amplitude can be established using the same methods as in the flat 



space case, as was shown in 12 



16 



18 



21 



22j . For pair production in an external elec- 



tromagnetic field, supposing that the fields are coupled by the elementary electric charge e, 
the expression of the transition amplitude is: 

A e - e+ = -ie I d A x[-g(x)} 1/2 U f , x (xtf ■ A(x)V f ,, x ,(x) (7) 



where e is the unit charge of the field and the fields Uf ) \(x) , Vp> \>(x), are supposed to be 



exact solutions of the free Dirac equation in the momentum basis 



4| . Using the fundamental 
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solutions of the Dirac equation given in equation (128]) from Appendix and the potential 
and setting Q = e, we finally arrive at the amplitude: 



/ppi 



16\p + p'\u 



+e" 



-sgn(XX')e- nk 
dzzHjS(pz)HjS(p'z) 



dzzEf\pz)Ef\p' z) 
Ca (#)(^ • n)rjx'(p '), 



(8) 



The result of the spatial integral {t| was included in the above amplitude. In the temporal 
integral a new variable z = —t c was introduced and sgn is the signum function. To solve these 
integrals we use the relation between Hankel functions and Macdonald functions jll[ 12 



, arriving in this way at the integrals discussed in Appendix. The result of our calculations 
can be expressed in terms of unit step functions 9, Euler gamma functions T and Gauss 
hypergeometric functions 2 -Pi: 



8n 2 m\p + p'\ 



-sgn(XX') (p- l e(p-p')f k 



+ p'- 1 6(p' -p)f k 



^(p)(a-n)r)x>(p'), 



(9) 



where m is the electron mass. The functions f k ^ ), which define the amplitude are given 
by: 




with the specification that the function f k yj^j is obtained when one makes the substitution 
p p', and the parameter k = m/oj is the particle mass/expansion factor ratio. When the 
momenta p, p' are equal, the result of our integrals can be obtained using the relations 
between Hankel and Bessel J functions and then the integrals can be solved in terms of the 
delta Dirac function 5(p — p'). As we show in |26( , the delta terms will cancel between them 
for each of the two integrals with Hankel functions. 

We will explore in the rest of the paper the physical consequences of our calculations. 
It is obvious from equation ([9]), that the only possible transitions are those in which the 
momenta of the electron and positron are not equal. From here we conclude that the law of 
conservation for the modulus of momentum is lost in this process. The background which 



enables the breaking of momentum conservation law is the external electric field, not the 



de Sitter geometry, because the geometry dD) has spatial translation invariance as an exact 
symmetry, so the momentum is conserved 26( . 

Further let us study the properties of our amplitude by drawing the graphs of the real 
and imaginary parts of fk (jf) (the analysis is similar for function fj. f^rj)> as function of 
parameter k = m/u, for different values of the ratio = x(— ~)- The parameter k encodes 
the influence of space expansion upon the pair production process. 




FIG. 1: The real part of fk(x) as a function of k. The solid line is for \ = 0-1 an d the dashed line 
for x = 0.9. 
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FIG. 2: The imaginary part of fk(x) as a function of k. The solid line is for x = 0.1 and the dashed 
line for \ = 0-9- 
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FIG. 3: The real part of fk(x) as a function of k. The solid line is for \ = 0.001 and the dashed 
line for x = 0.00001. 
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FIG. 4: The imaginary part of fk(x) as a function of k. The solid line is for \ = 0.001 and the 
dashed line for \ = 0.00001. 



From our graphs we observe that both the real part and the imaginary part of the function 
fk are nru te i n origin and very convergent for large values of k (see Figs. (pQ)- ([2])). As 

the ratio of the momenta takes small values, we observe that these functions become 
oscillatory (see Figs. (JSJ)- (SJ) ) - This oscillatory behavior is the result of the behavior of Gauss 
hypergeometric functions as their algebraic argument 1 — (jfrj approaches to one (or jj 
approaches zero), combined with the oscillatory factors yyj l ~ lk . Squaring the amplitude 
and summing after final helicities A, A', we obtain the probability of pair production. 

We find that the probability of fermion pair production in an electric field is: 



5 EM. 



IE 



2 ^— ^ 647r 4 m 2 \p + p 



- \Z+(p ) (a ■ n) Vy (p >)\ 2 {p >-H{p ' - p) 



+ p 9(p — p 



U[ p - 

p 



sgn{X\')[f k ll> 



*2 P 



}• 



(11) 



P J \p 

The total probability is obtained integrating equation fill I) after the final momenta p,p': 
P*°* e+ = J V e - e + (f^3-jf^3- From the general formula for probability (JTT1) . we observe that 
in the process of pair production by an electric field the helicity is conserved when A = — A' 
and the law of helicity conservation is broken when A = A'. We must specify that the 
non-conservation of helicity is driven by the fermion mass and not by the de Sitter geometry 
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26[ . Plotting our probability ( ITT]) as function of parameter k, we obtain the results from 
Figs. 
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FIG. 5: V e - e + as a function of k for x = 0.9. The dashed line represents the case of helicity 
conservation and the solid line represents the case when helicity is not conserved. 




FIG. 6: V e - e + as a function of k for \ = 0.1. The dashed line represents the case of helicity 
conservation and the solid line represents the case when helicity is not conserved. 
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FIG. 7: V e - e + as a function of k for x = 0.001. The dashed line represents the case of helicity 
conservation and the solid line represents the case when helicity is not conserved. 




FIG. 8: V e - e + as a function of k for x = 0.00001. The dashed line represents the case of helicity 
conservation and the solid line represents the case when helicity is not conserved. 



10 



From Figs.(JS])-(JHl) we observe that the electric field of the type considered here can produce 
pair of fermions only in the early universe, when the expansion factor was sensibly larger 
that the particle mass (oj > m). This result confirms the predictions from [s-lO] , were it 
was proven using a WKB approach, that the rate of pair production was important only 
in the early universe. As x — )• 0, the probability becomes oscillatory in both conserving 
and non-conserving helicity cases. It is also important to observe that the probability for 
a helicity non-conserving process seems to be sensibly larger that the probability for a 
helicity conserving process for % — > (see Fig. ((§])). Also from Figs. dSD-©, we observe 
that it is more likely to produce pairs of fermions, which have the ratio of the momenta 
X very small. To conclude, the presence of an electric field in this geometry will favour 
the processes in which for example ,the electron momenta will be small comparatively to 
the positron momenta, with the specification that large momenta don't refer here, to the 
relativistic case. As x is close to unity, we observe that the probability of pair production 
become sensibly smaller (see Fig. (jSJ)), than in the case of small x ( see Figs. (JZJ),©)- We 
also observe that our probability is proportional with a factor \p + p'\ 2 , from which it is 
clear that in the limit of large momenta our amplitude/probability will vanish. From here 
we conclude that only fermions with small momenta were produced in the presence of an 
electric field, due to the early expansion of the universe. 

Let us do a more detailed analysis in the helicity space. We use helicity bispinors 
£a(p) >^a'(p 0> for A = A' and A = —A', and restrict the analysis only to the square of 
our amplitude. In an orthogonal local frame {e*j}, we take the electron and positron mo- 
menta in the plane (1, 2), denoting their spherical coordinates as p = (p,a, (3) , p ' = (p 7, 5) 
and we consider that the unit vector which gives the direction and orientation of the electric 
field, has the spherical coordinates ft = (1, 9, ip) , where a, 7, 6 G (0, 7r); (3, 5, (p E (0, 2n). 

Corresponding to the two cases of conserving/nonconserving helicity in the pair produc- 
tion process by an electric field, the bispinor product £^~(p)(<? • n)r)\>(p') can be calculated. 
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In the case of helicity conservation this gives: 

£t{p){& ■ n)T]_i(p') = cos ^— j cos cos(6>) + e l ^~^ sin ^— J cos sin(6 l ) 
+e i(<5 "^ cos (|) sin sin(0) - e^"^ sin (|) sin cos(0), 
{p){B ■ n)r)i(p') = cos ^— j cos cos(6>) + e~ l ^ + ^ sin ^— j cos sin(6>) 
+e^" 5) cos (|) sin sin(0) - e- i(5+/3) sin (|) sin cos(0). 

(12) 

For the helicity nonconserving case the result is: 

ft(p)(5 ■ 7i)?7i(p ') = e" w cos sin cos(0) + e 4 ^"^ sin (|) sin sin(0) 

+e~^sm (|) cos (|) cos(0) - e-^cos (|) cos (|) sin(0), 
^x(p)(a ■ n)r)_±{p') = -e i<5 cos (|) sin Q) cos(0) - sin (|) sin sin(0) 

-e~ i/3 sin (|) cos cos(0) + e** cos (|) cos sin(0). 

(13) 

Taking now /3 = 7r,<5 = y? = and squaring the bispinor product we obtain: 

' cos(0) sin (2=2) + sin(0) cos (2=2) | 2 , for A = A'. 



(14) 



From equation ( Fl4|) if we set a = tt, 7 = = 0, corresponding to the situation when the 
electron-positron pair is generated on the direction of electric field, but they move in the 
opposite senses, we observe that the probability in the helicity conserving case is zero, while 
the probability of a process where helicity is not conserved becomes maxim. From here 
we can conclude that only the production processes that don't conserve the helicity could 
produce pairs that move in opposite senses, increasing in this way the chances of separation 
between matter and anti- matter. Now if we set a = 7 = 9 = 0, corresponding to the 
situation when the electron-positron pair is emitted on the direction of the electric field, but 
they move in the same sense, we observe that the probability in the helicity conserving case 
becomes maxim and is zero in the helicity non-conserving case. In this case, it is more likely 
that the pair will annihilate than become separate. Also from (I12 p .(ll3 p . p3)l . we observe that 
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the electron-positron pair could be emitted in other directions, which make various angles 
with the direction of the electric field, but with smaller probabilities. 

To conclude, in both helicity conserving/non-conserving cases, the most probable transi- 
tions are those where the pairs are emitted on the direction of the electric field. 

Let us further study the limit cases of our result. In the null mass case it is obvious that 



our amplitude/probabi 
other methods in 



ity vanishes. This is a confirmation of the results obtained using 



281 ] , where it was proven that there is no production of null mass 



fermions in the early universe because of the conformal invariance. 

For large k we can use r(|+i/c)r(| — ik) = (| — ik) ^^7^ and 2 Fi (a, b; b; z) — (1 — z)~ a 
to obtain an approximation for the function in the limit of large k. The final result is: 



/ fc »i( X ) ^ 2k (J - ik) ne^ k x^. (15) 

This function is very convergent and goes rapidly to zero for large k. This also can be seen 
from our graphs (see Figs. (CD)- (j4])). In the Minkowski limit k — > 00 (m >> u>), our analytical 
calculation gives = 0, confirming in this way the results from our graphs mentioned 
above. As k goes to large values, our amplitude vanishes. In Minkowski QED 

aai 



1 ik 



, pair production by the electric field of a charge, as considered here, is forbidden. This is 
because this field is too weak to produce pairs and the law of conservation for energy will 
be broken. The large mass approximation is also good in considering pair production for 
the present day expansion, this is because of the ratio mc 2 /hw ~ 5 ■ 10 37 , (where m is the 
electron mass). In the large mass case our probability is proportional with the factor e~ 2n ^ 
( see (fl5|) ). which makes the probability negligible for m » uj. The conclusion is that there 
is no fermion production in the electric field of a charge for the present day expansion , 
the probability being large only in the early universe, when the expansion factor was much 
larger than the particle mass. To conclude, the phenomenon of pair production, would be 
observable today only in strong electric fields. 



IV. NUMBER OF FERMIONS 



Let us translate our perturbational result in terms of density number of produced parti- 
cles. In the vast majority of the papers related to the particle production, the outcome of 
the calculations is the density number of created particles. So our reader could ask how we 
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will compute the density number, taking into account that we use a perturbative approach. 

First let us comment on the previous results related to the density number of created 
particles. The most important result in this direction was obtained by L.Parker using 
a WKB approach 10]. This result shows for the first time that the density number 
of created particles was important only in the early universe, which corresponds to the 
inflation phase. However there are results concerning the density number of created particles 



obtained using the method of Bogoliubov transformations. As was shown by Fulling |30|, 
there are many cases of interest where this method gives divergent results. This happens 
because the density number of created particles is not a function of the final momenta and 
as a result of that, the integration after the final momenta will give a quantity which is 



20|, 



linearly divergent. This important observation is valid for the results presented in 
where an integration after final momenta gives an infinite number of particles. So we can 
not even speak about a regularization of the integral simply because we do not have any 
dependence on final momenta 1_9|, |20] . In addition, when we compute the number of created 
particles in an expanding universe, this will be in general a time dependent quantity. This 
is the result of the fact the the volume will expand in time. For example, in de Sitter 
case discussed here, the physical spatial coordinates are x? h = a^e"*, and from here we 
observe that volume will expand as V P h = V e 3ujt . The above mentioned problems seem 



29] the density number of 



to receive little attention despite their crucial importance. In 
created particles was computed in a Robertson- Walker metric, and this is one of the few 
papers which takes into account the effect of the field interactions upon particle production 



processes. Moreover in 29], the correct density number of created particles was computed 
taking into account that the volume is a function of the expansion factor. 

Here we give the main steps for computing the number of created fermions in our case. 
Squaring the amplitude and summing after the final helicities A, A' we obtain the probability 
of pair production ( TTTT) . Then following the result from 29], the ratio probability/ volume, 
is the density number of fermions : 

A/-= — = -V |A ' e+|2 . (16) 
V ph 2^ Ve™ { ' 

This quantity can be computed now explicitly using the equation for probability pip . The 
total number of produced fermions can be obtained integrating the density number above, 
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after the final momenta p, p'\ 

Kftot _ [J^J^d^_ 

" V ph J V e 3 ^ (2vr) 3 (2tt) 3 ' 1 ' 

Here a few comments are in order. First our density number of produced fermions defined 
in (|T6|) . is dependent on the final momenta p, p' . However the integrals that result are 
not known in mathematics, these being a combination of two hypergeometric Gauss func- 
tions and a fraction of the type g^?p , where u, \i are complex numbers. Despite of these 
difficulties we will further give a way to calculate the total number of fermions. 

First from our graphs for probability Figs. flH])-©, we observe that there are larger 
probabilities to produce pairs of fermions which have the ratio of the momenta very small 
X << 1- This is the case when for example the positron momenta is large comparatively 
with the electron momenta (p' » p). Then for p' » p, we can rewrite equation (fTU]) . 
using (jSTj), as follows: 



fk [ — 



l-ik 




(18) 



Now we can simplify the hypergeometric functions, 2-^1(0, b;b; z) — (1 — z) a . Then using 
r(z)r(l — z) = s ^ z ^ and |T(l/2 + ik)\ 2 = , the final expression of the function f k is: 

P \ k TT 



p 



cosh(7r/c) - 



ifc / \ 1— ik 
V \ IP 



p' J \p' 



(19) 



The total probability will be calculated for the helicity nonconserving case and we will 
consider the most probable transitions which produce pairs that move on the direction of 
electric field. In the previous section we see that only the helicity nonconserving processes 
could produce pairs that are separable. This is equivalently to take pp' = \p\\p'\ cos(a + 
7) = Hp Hp 'I? corresponding to an angle between momenta vectors equal with tt (helicity 
nonconservation) and from here \p + p'\ 2 = p' 2 — 2pp' + p 2 . For A = A' the summation after 
helicity bispinors become: 

^t(p)(^-n)Vx(p')\ 2 = mp)(a-H)m(pT + \^(p)^-n)V-dp')\ 2 (20) 
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and if we consider that the most probable transitions in the helicity nonconserving case are 
obtained when a = n, 7 = = 0, corresponding to the situation when electron and positron 
move in opposite senses, the sum after helicity bispinors (l2Up . simplify to a numeric factor. 
For the momenta modulus integrals we consider that p' G (0, oo) and we restrict the limits 
of integration for p to p G (0, A). This cutoff for the electron momenta p is in accord with 
our suppositions that p « p' . 

Using (Tl9]) we compute the quantity that defines our probability: 



k 2 n 2 



pi 



+ 2-2 

pi 



-2 A 



1 - 

l+2ife 



cosh 2 (7rfc) 



2ik 



2-2ik 



P_ 
P' 



-2ik 



2+2ik 



(21) 



For p' » p the ratio ^ ~ 0, and only the term that don't contain the ratio p/p' and the 
two terms with purely imaginary powers from (1211) will give contribution, the result being: 



k 2 7l 2 



cosh (irk) 



pi 



2ik 



-2ik 



(22) 



Then the total probability for production of electron-positron pairs that move in opposite 
senses (helicity nonconserving case), but on the direction of the electric field vector, in the 
case p' » p is: 

e 4 fc 2 f dVL r 



V 



x 



tot 



32ir 2 m 2 cosh 2 (iik) J (2vr) 3 



dp' 



(p' 2 — 2pp' + p 2 ) 



2- A 



dVtpi 

2ik 



p 2 dp. 



(23) 



The integrals that help us to establish the expression for total probability are given in 
Appendix. Completing the momentum integrals calculation, the final result for the total 
probability turns out to be: 

a 2 k 2 A 2 



V 



tot 



■ cosh(27rA:) _ i 
87r 4 m 2 cosh 2 (vrA;) L ; sinh(27rfc) 



(24) 
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where a = e 2 /(47r) is the fine-structure constant. Here a few comments need to be made. 
We know that the electric field of the type considered here doesn't produce pairs of fermions 
in Minkowski QED. Our result for total probability in equation (I24p . vanishes in the flat 
space limit (k = oo). To see that, we plot the functions dependent of k from equation 



g(k) 



cosh (irk) 



(2vrfc) 



cosh(27rfe) 
sinh(27rfe) 



0.1-1 



0.06- 



"oh 



0.04- 




FIG. 9: g(k) as function of parameter k 



The conclusion is that the total probability for pair production in an electric field of a 
charge, is nonvanishing only in the early universe as we predict with our graphs for proba- 
bility in Figs. dSJ)- (JHI) - Then the final result for total number of fermions in volume unit will 
be: 



a 2 k 2 A 2 



tot 



87r 4 m 2 cosh (irk) 



, , cosh(27rfc) 



-'Scut 



(25) 



sinh(27r&;) 

We see that the total number of fermions is a time dependent quantity and decreases expo- 
nentially in time. This is natural since the expansion factor decreases also with time. From 
here we can conclude that the effect presented in this paper could be observable only in 
strong gravitational fields. Another interesting observation is that the total probability and 
total number of fermions depend on the parameter A, which establishes the upper limit of 
the electron momenta. 
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V. CONCLUSIONS 



We compute the amplitude/probability for pair production by the electric field of a charge 
in the de Sitter expanding universe. The probability of pair production is nonvanishing only 
in the early universe when the expansion factor was large. From our calculations we recover 
the Minkowski limit, where the amplitude and probability for this process vanish. We found 
that that the most probable transitions are those in which the pairs are emitted on the 
direction of the electric field. 

We also translate our perturbational calculations in terms of density number of fermions. 
A method to compute the total number of fermions in the case when the ratio jj is close 
to zero, was proposed. The total number of produced fermions was obtained and we prove 
that is nonvanishing only in the early universe. 



VI. APPENDIX 

Let us introduce unit normalized helicity spinors [6 ; ] for an arbitrary momentum vector 

V ■ £a(p) and Vx(p) = ic2[£\(p)]* 

apUP) = 2pXUp) (26) 

with A = ±1/2 and where a are the Pauli matrices and p =\ p\. The particle spinors have 
the form 



-pi+ip2 



^W^(^J< 7 1 (27) 

Then, the positive/negative frequency modes of momentum^ and helicity A derived in 
4( , assuming gamma matrices in Dirac representation, are [l, % 5] : 

where H^\z) are Hankel functions of the first kind and k = —,u± = \ ± ik. The 
negative frequency modes will be obtained via charge conjugation operation Vp t \(t,x ) = 
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For solvin g ou r integrals with Hankel functions we use their relations with Macdonald 



functions llll, 112L [23j : 



H^ 2 \z) = T (^j e^' 2 K v { T iz). (29) 



In this way we arrive at the integrals of the type 23J : 



dzz A K^{az)K y {bz) = — r 

1 [i - A) 



'l -X + fj,-u \ / l-X-fj,-u \ f l-X + fi + v 1- \- [i + v b 2 
2 J V 2 y 21 V 2 ' 2 ' ' a 2 

Re(a + b) > , Re(\) < 1 - \Re(fj) \ - \Re{v)\. ( 



In our case A = — 1 and the second condition for convergence is satisfied. We also observe 
that in our case a, b are complex and for solving our integrals we add to a a small real part 
a — > a + e, with e > and in the end we take the limit e — > 0. This assure the convergence 
of our integral and will correctly define the unit step functions and /& functions. 

For the computation of the total number of particles we use the next relation between 
hypergeometric functions: 

2 F 1 (a, b;c;z) = — — — 2 Fi(a, b ; a + b - c + 1; 1 - z) 

l (c — a)l (c — o) 

+ (1 - zf- a - b T(C) 5\iL~ C) 2F 1 (c-a,c-b;c-a-b+l;l-z). (31) 

r(a)r(6) 

The integrals that helps us to calculate the total probability are: 

=F2ifc 

, (32) 



dp' 1 f°° dp'p' ±2ik - 2 ^ k (- l P 



(p' 2 — 2pp' + p 2 ) p 1 J (p' 2 — 2pp' + p 2 ) p sinh(7rfc) 
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